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ABSTRACT 

A  new  semi -empirical  expression  has  been  developed  which  predicts 
the  evaporation  rates  of  pure  liquid  hydrocarbons  under  a  wide  range  of 
conditions.  A  code  has  been  formulated  which,  using  the  above  results, 
calculates  the  evaporation  rates  of  binary  mixtures  of  hydrocarbons. 
Extension  of  these  results  to  real  fuels  is  possible  by  finding  the 
appropriate  binary  mixture  which  approximates  the  fuel .  Analyses  of 
the  flammable  envelopes  surrounding  pools  of  pure  hydrocarbons  and 
fuels  have  also  been  performed. 
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A  wide  range  of  vaporization  problems  both  experimental  and  theo¬ 
retical  have  been  explored.  Sections  I  -  V  present  theoretical  analyses 
of  diffusion  vaporization,  v;  -ization  under  uniform  flow  conditions, 
and  vacuum  vaporization.  These  sections  provide  the  basis  for  the 
semi -empirical  theoretical  models  of  the  more  practical  problems. 
Practical  experimental  results  as  well  as  scmi-empi rical  theoretical 
analyses  cc  these  results  are  presented  in  Sections  VI  -  X.  These 
sections  should  be  of  greatest  interest  to  persons  engaged  ‘n  field 
tests . 

Vaporization  even  under  pseudo-steady  state  conditions  is  a  compli¬ 
cated  p'  enomena.  Vaporization  rates  for  pure  liquids  are  dependent  on 
the  air  velocity  distributions  above  the  liquid,  liquid  surface  tempera¬ 
ture,  vapor  pressure,  diffusion  coefficient,  molecular  weight,  and 
geometric  factors.  A  semi -empirical  expression  has  been  developed 
(Section  VI)  which  includes  these  parameters  and  which  accurately  models 
evaporation  from  an  open  pool  into  a  laboratory  room  where  there  is  slow 
collective  air  motion  (average  speed  -  3  cm/sec)  with  velocity  fluctua¬ 
tions  (see  Figure  l) .  This  velocity  distribution  is  closely  reproduc¬ 
ible  over  the  period  of  the  evaporation  rate  measurements .  Time  averaged 
evaporation  rates  measured  under  these  conditions  are  reproducible  to 
better  than  ±  10%,  With  slight  modification  the  expression  should  pre¬ 
dict  vaporization  rates  under  a  wide  range  of  flow  conditions  not  just 
those  mentioned  above. 

The  evaporation  rates  of  simple  mixtures  are  dependent  on  all  of 
the  above  pure  liquid  factors  and,  in  addition,  to  the  liquid  phase 
mole  fractions  and  to  the  activity  coefficients.  A  simple  binary 
mixture  (1/3  pentane  -  2/3  octane  by  volume)  has  been  discovered  which 
simulates  the  vaporization  characteristics  of  a  batch  of  combat  gaso¬ 
line.  A  somi-empirical  theory  (Section  VII)  which  predicts  the  vapori¬ 
zation  rate-time  curve  for  a  binary  mixture  was  developed  and  coded. 

The  code  predicts  the  vaporization  rate  change  with  time  very  accurately 
for  a  i/3  pentane  -  2/3  octane  mixture  (under  the  same  laboratory  room 
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flow  conditions  described  for  pure  liquids) . 

In  conjunction  with  evaporation  rate  analysis  the  associated  problem 
of  concentration  profile  development  around  liquid  pools  was  explored 
both  for  pure  liquids  and  mixtures  (Section  VIII).  Profiles  are 
determined  with  a  spark  probe  by  locating  the  flammable  volume.  Data 
on  steady  state  and  time  development  of  flammable  volumes  are  presented. 
A  limited  theory  which  predicts  the  height  of  the  flammable  volume  is 
suggested. 

A  graphical  analysis  of  the  heat  transfer  problem  associated  with 
an  evaporating  pure  liquid  is  developed  (Section  IX) .  This  analysis 
demonstrates  in  principle  how  the  liquid  evaporation  rate  may  be  de¬ 
termined  if  the  surface  temperature  is  unknown.  If  the  surface  tempera¬ 
ture  is  unknown  then  the  thermal  properties  of  the  air,  liquid,  and  base 
material  as  weiJ  as  the  heat  of  vaporization  must  be  known  in  addition 
to  all  the  "ther  properties  previously  mentioned  for  pure  liquids. 

The  discussion  thus  far  has  been  limited  to  pseudo-steady  state 
conditions.  Some  one-dimens io.isl  time  dependent  pure  diffusion  systems 
are  discussed.  Also  steady  state  vaporization  problems  in  one,  two 
(with  flow) ,  and  three -dimens ions  are  discussed.  These  problems  pro¬ 
vide  a  basis  for  theoretical  analysis  of  the  more  practical  problems. 

Although  somewhat  removed  from  the  rest  of  the  report,  a  section  on 
vacuum  vaporization  is  included.  There  is  much  literature  on  this 
subject  and  it  was  thought  that  a  section  devoted  to  this  topic  would 
help  put  it  in  proper  perspective. 

II.  VACUUM  VAPORIZATION 

Evaporation  rates  in  low  pressure  vaporization  are  often  less  than 
predicted  from  kinetic  theory  by  a  factor  a,  the  evaporation  coefficient. 
Much  attention  has  been  given  to  determining  a  and  numerous  methods  have 
been  employed. 1  * 2» 3*  These  studies  are  useful  in  the  design  of  molecu¬ 
lar  stills. 

References  are  listed  on  page  62. 
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Let  us  consider  vaporization  into  a  nearly  evacuated  container  of 
volume  Vc  from  a  liquid  of  surface  area  A^.  The  time  development  of 
the  number  density  of  molecules  (n)  in  this  situation  can  be  written 

AL 

n  =  —  a  c'  (ng  -  n)  (1) 

c 

where 

a  is  the  evaporation  coefficient 
c"  is  the  mean  velocity  of  the  gas  molecules  » 

Ts  is  the  liquid  surface  temperature 

M  is  the  liquid  molecular  weight 

R  is  the  universal  gas  constant 

ng  is  the  equilibrium  number  density  (number  of  molecules 
s  per  unit  volume  at  equilibrium) . 

The  above  equation  is  valid  only  where  temperature  changes  in  the  liquid 
can  be  ignored,  a  fact  which  restricts  its  application  to  liquids  of 
very  low  vapor  pressure  typically  glycerol.  Integrating  the  above 
equation  results  in 


n  -  n^  A. 

,  s  o  L  .  _ 


(2) 


where 

nQ  is  the  number  density  at  time  zero  t  *  0 

t^  is  the  time  required  for  the  number  density  to  reach  n^. 

Now  let  us  cover  the  exposed  liquid  surface  with  a  diaphram  which 
is  pierced  by  a  tiny  hole  of  area  A  .  The  number  density  in  the  space 
between  the  liquid  surface  and  the  diaphram  is  designated  n'.  The 
following  equation  can  be  written,  similar  to  Equation  (1)  for 
vaporization  into  Vc 


n  *  —  c'  Cn • 


n) 


(3) 
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integrating  this  equation  results  in 


fi-n 


n'  -  n 
_ o 

n'  -  n^ 


C4) 


A  fairly  good  assumption  is  that  n'  =  ng  for  very  small  holes.  The 
time  tj  is  the  time  required  for  the  number  density  to  reach  n^ ,  much 
longer,  of  course,  than  t^  above.  Thus  taking  the  ratio  of  Equation  (2) 
and  Equation  (4)  the  evaporation  coefficient  is  obtained 


o 


A  t ' 

A0  h 


(5) 


The  basic  method  described  above  and  known  as  the  molecular  effusion 
method  has  been  refined  by  Heideger  and  Boudart2,  The  best  value  of  a 
for  glycerol  appears  to  be  about  0.05. 

It  is  emphasized  here  that  a  appears  to  be  important  only  in  low 
pressure  vaporization.  It  represents  a  kind  of  surface  resistance  to 
evaporation. 

The  maximum  rate  of  vaporization  from  kinetic  theory  is  given  by 


(n) 


max  th. 


(6) 


The  experimental  maximum  rate  is  modified  by  a  and  is 


c'  A, 


(n)  =  a 

^  'max  exp. . 


—  ns  ’ 


(7) 


It  is  interesting  to  note  that  Equation  (1)  can  be  rewritten 
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n  -  n 


n  =  a  (L  c*)  (— — x— )  (8) 

-  L 

c 

where  L  =  V  /A.  represents  a  mean  free  path  length . 

C  C  L 

a  •  L  •  c'  has  the  dimensions  of  a  diffusion  coefficient.  Equation 
(8)  then  closely  resembles  the  diffusion  equation 


Vs  D 


but  in  Equation  (8)  the  effective  D  (a  •  L  ' '•> c')  is  very  very  large  due 
to  the  very  large  mean  free  path  length.  Cp  in  Equation  (9)  is  the 
vapor  concentration.  Thus  if  a  liquid  is  allowed  to  evaporate  into  a 
vacuum,  equilibrium  will  be  attained  much  faster  than  if  it  is  allowed 
to  evaporate  into  air. 


HI.  ONE  DIMENSIONAL  DIFFUSION  VAPORIZATION 
A.  Steady  State 

The  evaporation  mole  flux  (Npz)  for  the  one  dimensional  steady  state 
problem  with  a  specified  mole  fraction  (XpQ)  at  the  liquid  surface  (at 
height  zQ)  and  another  specified  mole  fraction  Xp^  (at  height  Zj)  at  the 
top  of  the  tube  is 


h  -  V 


1  -  X- 


where 

c  is  the  total  molar  concentration 
D  is  the  binary  diffusion  coefficient  of  F  in  A, 

Using  vapor  pressure^  instead  of  mole  fractions,  the  mass  flux  per 
unit  area  is  found  to  be 
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(11) 


The  above,  theoretical  expression  has  been  compared  to  experiment  for 
ID  open  tubes  (exposed  to  a  laboratory  room)  for  both  cyclohexane  and 
ethyl  alcohol .  .In.both  cases  it  is  assum'd  that  Ppi  is  zero,  that  is, 
that  the  air  motion  in  the  room  is  sufficient  to  maintain  a  nearly  zero 
concentration  at  the  top  of  the  tube.  In  addition  to  comparing  the 
above  one  dimensional  theory  with  experiment  the  results  are  compared 
with  the  following  mass  evaporation  rate  expression  inferred  from 
diffusion  vaporization  and  from  2D  boundary  layer  theory  (flow  theory) . 
This  expression  will  be  discussed  in  more  detail  later  in  this  report. 


m  = 
v 


=  4  (1  ♦  reCRo)1^2  (Sc)1/3) 

^  f  i  .  71  172  173 


1  +  J.  (i  +  2i*8(Re) 


(Sc) 


in 


where 

m 

V 

K* 

8 

Re 


*  total  mass  evaporation  rate 

«■  1.49 

*  0.34 

2  u„  a 
o 

s  .  ——-.i. .  - 

V 


RT 


•  D  •  a  (12) 


SC  a  V/D 

u  *  average  convective  air  speed  in  the  room  (taken  to  be 

0  3  cm/sec  in  this  case) 

v  #  kinematic  viscosity  (u/p) 

a  *  pan  radius 

D  ■  binary  diffusion  coefficient 
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PFo  =  vaPor  Pressure  at  liquid  surface 

=  atmospheric  pressure 

T  =  absolute  temperature  (°K) 

Mp  =  molecular  weight 

i  =  depth  of  liquid  below  pan  rim. 

The  experimental  results  obtained  here  at  BRL  and  the  results  from 
the  two  theories  are  compared  in  Table  I.  All  results  are  for  a  =  3.7  cm. 


Tabx.  . 

Finite  Open  Tube  Evaporation  Rates 
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T 
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ID.  THEORY 
RATE 

(«4 

min' 

FLOW 

THEORY  RATE 

CYCLOHEXANE 

0 

294°K 

.25 

co 

.23 

1 

297°K 

.023 

,024 

.022 

3 

297°K 

.0086 

.0080 

.0077 

;.3 

297°K 

.0053 

.0045 

.0045 

ETHANOL 

0 

295®K 

.085 

CO 

.096 

l 

296  ®K 

.0110 

.0102 

.0093 

2 

296  #K 

.0060 

.0051 

.0049 

3 

296°K 

.0039 

.0034 

.0033 

One -dimensional  open  tubes  (or  Stefan  tubes)  aro  one  of  tho  primary 
sources  of  gas  phase  diffusion  coefficients.  Mass  evaporation  rates 
from  the  tube  are  measured  and  the  one -dimensional  theory  is  used  to 
calculate  the  diffusion  coefficient  (D) .  A  carefully  controlled  flow 
of  air  (or  other  carrier  gas)  is  passed  over  the  mouth  of  the  tubes  to 
insure  zero  concentration  at  the  top.  Such  methods  have  been  used 
rocently  by  lee  and  Wilke4  (1954)  and  by  Altshullor  and  Cohen5  (1960) 
to  determine  diffusion  coefficients.  The  concentration  profiles  in  a 
Stefan  tube  have  been  investigated  experimentally  by  lleinzelmann, 

Wason,  and  Wilke6  (1965)  and  have  boon  found  to  agree  well  with  the  theory. 
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B .  Time  Dependent  Infinite  Open  Tube 

Exact  solutions  exist  for  the  one-dimensional  time  dependent  system 
of  infinite  length.  Discussion  of  this  specific  evaporation  problem  is 
presented  by  Bird,  Stewart  and  Lightfoot7  and  for  a  general  class  of 
problems  of  this  type  by  Danckwerts8. 

In  the  discussion  in  Bird  the  liquid  surface  is  kept  at  a  fixed 
level  (z  *  0)  and  it  is  assumed  that  the  flux  of  air  at  this  point  is 
zero  NAq  =0.  The  following  species  equation  is  solved 


3xp 

W 


32 


z 


(13J 


subject  to  the  following  boundary  and  initial  conditions 


xp  =  0  t  -  0 

XF  ’  xFo 

Xp  »  0  z  , 

The  above  equation  is  converted  to  dimensionless  form  and  the  solution 
is 


where 


v  _  1  -  erf  (Z  - 

X  *  rV  erf  V' 


(14) 


X  ■  *p/*p0 


Z  ■ 


/TuT 


1  Fo  dX 
21  "  xFo  ^ 


z  •  0. 


25 


<j>  may  be  thought  of  as  a  dimensionless  bulk  velocity,  in  fact  the  bulk 
molar  velocity  (v*)  is  related  to  $  by 

v*  »  <j>  /o7t 

v*  is  zero  only  for  xpo  =0.  v*  is  known  as  the  Stefan  flow9.  There 
can  be  no  diffusion  vaporization  without  some  convection  v*.  The 
solution  X  for  various  values  of  XpQ  is  shown  in  Figure  2. 

The  rate  of  evaporation  (volume  rate)  is  given  by 

d  VP 

*  A  /07t  (16) 

whore  A  is  the  cross  sectional  area  of  the  tube.  Relations  of  the 
above  type  have  boon  known  for  many  years,  in  fact  a  verification  of  it 
was  presented  by  M.  P.  Vaillant  in  191010  (Vaillant  uses  the  Pick’s 

Xp 

second  law  approximation!  4  small  and  $  %  — -  in  equation  (16))* 

/r 

It  is  of  interest  to  note  that  the  time  dependent  finite  open  tube 
problem  with  zero  concentration  at  the  top  and  fixed  concentration  at 
the  liquid  surface  (initial  assumption  of  zero  vapor  concentration  in 
the  tube)  has  an  analytical  solution  in  terms  of  Fourier  series.  The 
results  for  the  analogous  heat  flow  problem  are  presented  by  Tarifa11. 

C.  One -Dimensional  Closed  System 

In  connection  with  the  problem  of  concentration  profiles  in  a  fuel 
tank  a  one -dimensional  time  dependent  analysis  of  a  closed  system  was 
performed.  The  initial  conditions  are  that  no  fuel  vapor  is  initially 
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present  in  the  region  above  the  fuel  surface.  Boundary  conditions  are 

that  a  fixed  concentration  exists  at  the  liquid  surface  (c_  )  and  that 

8  c  F0 

the  flux  of  fuel  at  the  top  of  the  tank  is  zero,  D  (-5 — •)  =  0. 

z  =  L 

c 

The  equation  to  be  solved  is 


32  c, 


(17) 


subject  to  the  above  mentioned  conditions.  The  solution  to  this 
equation  is  obtained  by  a  standard  separation  of  variables  method  and 
a  Fourier  series  expansion  of  the  concentration.  The  result  is 


cp  ( z ,  t) 


CF< 


1  - 

“  nit 

n  «  1,3,S,7 
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rntz 

2L 


2  2 
n  it 
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It  is  important  to  note  that  the  above  result  is  only  valid  for  small 
fuel  concentration  at  the  surface  (certainly  no  greater  than  twenty  per¬ 
cent  by  mole  fraction). 

An  approximation  to  the  pressure  as  a  function  of  time  in  this 
closed  system  can  be  obtained  by  integrating  the  above  expression  for 
the  concentration.  The  result  is 


p*  Wo 


V ft 

n  «  1. 3,5,7 


2  2 

n  a 
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(id) 


where 


PA  is  atmospheric  pressure 

Ppo  is  the  vapor  pressure  of  the  hydrocarbon  fuel . 

The  concentration  profiles  as  a  function  of  time  and  the  pressure-time 
relationship  are  shown  in  Figures  3  and  4. 


IV.  TWO  DIMENSIONAL  VAPORIZATION  INTO  FLOW 

The  following  equations  are  considered  in  this  analysis  of  the  flow 
problem  whore  the  liquid  surface  is  considered  flat  and  rigid, 

(Figure  5) . 
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with  the  following  boundary  conditions 
for 

*  «  0  or  y  *  T  *  T*. 

at 
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figure  3.  Contemrcrtion  Profiles  «n  Closed 
Fuel  Tank . 


No  chemical  reactions  are  assumed  here  and  species  A  represents  ambient 
air  into  which  species  F  is  evaporating.  The  following  properties  are 
assumed  constant  p,  u,  C^,  k,  c  and  D.  First  the  continuity  equation 
is  solved  for  u^. 
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dy 
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or  in  terms  of  the  molar  flux  of  F  at  the  surface  (NpQ) 


(22) 


Fick's  first  law  gives 
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3xp 
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(23) 


y=0 


but  NAq  »  0,  thus  substituting  for  Npo  into  equation  (22)  above  the 
result  is 
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y=0 


but  the  original  assumption  p  =  const.,  c  *  const.,  force  MA . *  Mp  and 


M  •  c 

thus  -  »  1.  Therefore, 
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y=0 


u  is  substituted  into  the  momentum,  energy,  and  species  equations  to 
obtain  the  following  common  form  for  the  three  equations 
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'  ^  *“  t0)  '  ^  2*u  in\  *'  = 


or  (defining  K) 


(26) 


-A  |  -  K  ♦  /£  2ttu  dn  |  tt  '  =  it" 


where  it  may  be 
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T  -  T 


TT_  = 
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and  A  may  be 


.  V  .  .  V  .  V 

Au  '  v  “  1‘  AT  ■  a  AFA  *  D 


and  the  prime  denotes  differentiation  with  respect  to 


=  Y  jS 

2  1  vx  * 


The  boundary  conditions  are 

n  =  0,  n  =  0;  n  +  ®»  ir  ■*  1 . 

Analytical  solutions  are  obtainable  in  the  2D  case  only  for  small  K  and 
moderate  to  large  values  of  Ap^. 

An  interesting  point  can  be  made  at  this  juncture.  Consider  the 
general  equation  written  for  the  momentum  case 

-  Au  K  *  2  'll  *>)  *i  =  *i’ •  (28) 

Now  with  K  »  0  (no  evaporation  or  condensation)  the  above  becomes 

-  {/q  2  \  “  *u'>  Au  “  l'  (29) 
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Now  letting  ir  * 


u 


~  the  above  is 
3ri 


-  2f 


32f 

3n2 
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(30) 


which  is  just  the  Blasius  equation12  and  has  a  well  known  solution. 

To  obtain  the  evaporation  mass  flux  the  species  version  of  equation 
(26)  above  is  solved  giving  (for  small  K) 


^FA  AFA*  ^ 


-  r  Q/3.  u) 

“  rTi/3,  “) 
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where 


1.328 


The  derivative  at  the  liquid  surface  is  (for  small  K) 


(32) 


*FA  C°»  AFA*  K)  =  °-68  AFA1/3*  (33) 

Now  using  the  definition  of  K  and  assuming  Xpm  =  0  the  result  is 
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Substituting  for  (0)  and  assuming  XpQ  small,  the  above  becomes 
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0.68. 


(35) 
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Pick's  first  law  gives 


Hp^  is  the  mass  flux  per  unit  area  of  fuel  at  the  surface.  Solving  for 
nc  ,  the  following  interesting  result  is  obtained 


The  above  result  was  obtained  assuming  Mp  =  M^.  Experiment,  however, 
indicates  it  is  valid  even  for  Mp  f  M^.  Including  the  molecular  weight 
factor  (Mp/M.)  the  result  is 


where  PFo/PA  “  xFo  and  PpQ  is  the  fuel  (liquid)  vapor  pressure  at  the 
surface,  is  atmospheric  pressure. 

The  above  expression  for  npo  can  be  integrated  over  a  surface  area 
to  give  the  total  mass  evaporated  per  unit  time  for  the  area  into  a 
unidirectional  flow  parallel  to  the  surface. 

To  integrate  npQ  over  a  circular  liquid  surface  the  surface  is 
considered  to  be  made  up  of  a  series  of  rectangles  with  long  sides 
parallel  to  the  flow.  The  two  dimensional  model  above  is  applied  to 
each  rectangle  assuming  no  interaction  between  rectangles.  The  follow¬ 
ing  diagram, illustrates  the  method  of  integration. 


The  total  mass  evaporated  per  unit  time  from  a  rectangle  of  width  'dy'  is 


Integrating  over  *x'  results  in 

n  2/3  PPn  7  7  l/4Mp 

mv  (dy)  -  0.34  (J-)  •  ♦  (uw  •  u  ♦  p)1/<4  2  /I  (a'V)  g-  dy.  (40) 

A  A 

The  total  rate  is  obtained  by  integrating  over  y  from  -a  to  +a.  The 
integral  is  obtained  using  Simpson's  rule  on  a  9810A  H.P.  calculator. 

The  result  is 

.  n  2/3  Ppn  Mp  i  m  x/o 

mv  "  °*34  ($  •  (u.  •  w  •  P)A/  •  4  /2  •  0.863  •  aA/*  (41) 

A  A 

or  written  in  terms  of  the  Reynold's  number  and  the  Schmidt  number 

P  NL 

my  »  0.34  *  Z2-  *  jp  ‘  (Sc)l/3  •  (Re)1/2  •  p  •  a  •  3.45  •  D  (42) 
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for  evaporation  from  a  circular  fluid  surface  into  a  uniform  stream  of 
air  parallel  to  the  surface,  where 

2  u  a 

\)  00 

Sc  =  jr  Re  =  - 

D  v 

p  is  the  density  of  air  and  vapor  nearly  equal  to  (  air  density) 
for  small  vapor  pressures.  A  square  with  side  length  '2a'  would  give 
a  closely  similar  result,  3.45  would  be  replaced  by  4.00. 

A  result  showing  a  like  dependence  on  the  Schmidt  and  Reynolds 
numbers  for  evaporation  from  a  planar  surface  has  been  obtained  by  other 
investigators13. 


V.  THRtL  DIMENSIONAL  DIFFUSION  VAPORIZATION 


One  problem  of  interest  is  evaporation  from  an  open  pan  into  a 
perfectly  still  atmosphere.  In  such  a  situation  diffusion  of  vapor  into 
the  air  controls  the  evaporation.  In  its  simplest  form  (assumption  that 
the  vapor  pressure  is  small  compared  to  atmospheric  pressure)  the 
problem  reduces  to  the  solution  of  Laplace's  equation  for  the  concentra¬ 
tion  in  three  dimensions  subject  to  the  following  boundary  conditions: 
zero  concentration  at  infinity,  tero  vapor  flux  from  r  *  a  (pan  radius) 
to  infinity  in  the  plane  of  the  pan,  and  equilibrium  vapor  concentration 
at  the  liquid  surface.  These  boundary  conditions  are  summarized  in 
Figure  6.  The  solution  curves  of  equal  concentration  are  ollipsoids 
centered  on  the  pan  and  are  shown  in  Figure  7.  The  solution  along  the 
z  axis  has  a  simple  form  and  is  given  by 
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Likewise  the  solution  in  the  plane  of  the  pan  has  a  simple  form 


CF 


2  c, 


Fo 


sin 


(44) 


36 


0 


'‘“LIQUID  SURFACE  IN 
*-y  PLANE 

0  S  r  S  a  Cf  »  CF< 


V2Cf*  0 

VAPOR  CONCENTRATION 


Figure  6.  Boundary  Conditions,  Three-Dimensional 
Laplace  Equation. 
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Figure  7  Solution  to  the  Three-Dimensional 
Laplace  Equation. 


For  moderate  to  large  values  of  r/a  the  equi -concentration  profiles  are 
very  nearly  spherical .  The  mass  evaporation  rate  from  a  pan  of  radius 
•a'  according  to  this  theory  is 
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where  all  the  symbols  have  been  defined  previously. 

The  above  solution  is  based  on  the  assumption  of  small  surface 
concentration  (small  PpQ/P^) .  The  evaporation  rate  expression  extended 


to  large  values  of  PpQ/PA  is 
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The  above  result  is  due  originally  to  Stefan  as  reported  by  Thomas  and 
Ferguson14  and  of  course  reduces  to  the  previous  result  for  small 

PFo/PA* 


VI.  COMBINED  MASS  EVAPORATION  RATE  EXPRESSION 

The  combined  mass  evaporation  rate  expression  was  formulated  to  fit 
the  data  of  a  laboratory  experiment  which  simulates  real  world  conditions. 
The  experiment  consists  of  a  liquid  in  a  circular  pan  which  is  allowed 
to  evaporate  into  a  room  where  there  is  a  slow  nit  flow  of  some  3  cm/sec. 
passing  over  the  pan  (average  velocity  component  in  the  direction  of 
flow).  The  fluctuations  in  the  air  flow  both  parallel  and  transverse  to 
the  flow  direction  are  some  t  S  cm/sec.  These  fluctuations  lead,  of 
course,  to  reversals  in  the  flow  direction  at  times.  DISA  (DISA  Elec¬ 
tronics,  Franklin  takes,  N.J.)  and  TSI  (The  mo -Systems  Inc.,  Saint  Paul, 
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Minn.)  low  velocity  anemometers  were  employed  for  the  flow  measurements. 
The  flow  velocity  distribution  (see  Figure  1)  is  quite  reproducible  over 
the  period  of  the  evaporation  rate  measurements.  The  experimental 
evaporation  rates  (time  averaged  over  a  few  minutes)  are  reproducible 
to  within  ±  10%  (nearly  steady  state  thermal  conditions  must  be  attained 
indicated  ty  nearly  steady  liquid  surface  temperature  before  measurements 
are  made) . 

The  theoretical  results  of  the  previous  two  sections  (diffusion  into 
air  flow  and  pure  diffusion)  can  be  combined  together  into  a  single 
equation  for  the  mass  evaporation  rate.  The  result  is  (assuming  a 
unidirectional  flow  over  the  pan  and  small  surface  concentration  of 
fuel  vapor) 

p  p 

my  -  4  (l  ♦  0.863  &  (Re)1/2  (Sc)l/5)  (^)  *  •  Np  •  D  •  •  .  (47) 

A 

An  equation  of  similar  form  showing  the  same  dependence  on  Reynolds 
number  and  Schmidt  number  (due  to  Frossling  as  reported  by  Fuchs9)  has 
been  used  for  years  to  successfully  describe  droplet  evaporation  even 
at  very  low  ReynoMs  number  (*  10)  where  there  is  no  current  theory. 

It  is  very  important  to  note  that  the  component  of  vaporisation  in  the 
above  equation  due  to  flow  is  greater  than  the  pure  diffusion  component 
even  at  flows  of  1  cm/sec  (for  the  range  of  pan  radii  considered) . 

Thus,  even  very  small  flows  significantly  effect  the  evaporation  rate. 

The  above  theoretical  Frossling  equation  is  valid  for  unidirectional 
flow.  The  actual  experiment  which  is  to  be  described  does  have  a  sisal  1 
unidirectional  flow  component,  however,  there  are  large  fluctuations 
about  this  velocity.  It  would  appear  to  be  important  to  obtain  velocity 
probability  distributions  along  three  mutually  orthogonal  axes.  These 
distributions  should  completely  characterise  the  flow.  External  flow 
is  very  important  in  the  actual  experiment.  The  flow  characteristics 
and  hence  the  evaporation  rate  cm.  be  changed  by  opening  and  closing 
air  conditioning  ducts  in  the  laboratory  room. 
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In  an  attempt  to  simply  describe  evaporation  into  the  laboratory 
room  (without  using  probability  distributions)  a  Frossling  like  equation 
was  adopted.  The  mean  measured  air  speed  (direction  not  considered) 
ove*-  the  liquid  surface  was  used  in  the  Reynolds  number  (-3  cm/sec) . 

The  equation  has  the  form 


»v  •  4  {1  ♦  V  8  (Re)1/2  (Sc)i/3)  In  ( 
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where  K*  is  a  constant  dependent  on  the  pan  geometry  and  the  flow 
velocity  distribution  (K*  is  determined  by  fitting  the  data  for  one 
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liquid,  in  this  case  cyclohexane),  a-  is  replaced  by  tn  C — y”! 

,  r© 

*  -  b — 

A 

in  the  hope  of  extending  agreement  to  higher  vapor  pressures. 

The  worst  deviation  between  theory  and  experiment  is  30\  with  other 
rates  in  much  better  agreement.  The  expression  fits  the  data  well  (see 
Table  II  to  compare  theory  end  experiment). 


The  question  of  the  effect  of  relative  humidity  an  evaporation  rates 
has  not  been  addressed.  Relative  humidity  does  not  appear  to  he  signi¬ 
ficant  for  ail  the  liquids  tested  except,  of  course,  for  water.  The  range 
of  humidity  encountered  in  these  experiments  is  from  4S  to  dd  per  cent. 

for  t/a  /  0  (t  is  the  depth  of  the  liquid  below  the  pan  rim)  t hr 
frossling  equation  above  has  been  further  modified.  The  modified 
equation  is  as  follows 
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for  very  large  t/a  the  above  reduces  to  the  one -dimensional  theoretical 
result  for  evaporation  from  am  open  tube  with  sere  intent  ration  at  the 


41 


top.  The  flow  in  the  room  maintains  approximately  zero  concentration  at 
the  end  of  the  tube  for  large  1/ a.  The  above  equation,  as  mentioned  in 
Section  IIIA,  was  applied  to  cyclohexane  and  ethanol  for  £/a  /  0. 

(See  also  Figure  8) . 

The  combined  rate  expression  has  been  shown  to  be  effective  in 
analyzing  the  experiments  conducted  thus  far.  The  expression  contains 
terms  which  are  characteristic  of  mass  transfer  in  flow  systems 

2/2  ]  /i 

(Reynolds  number  and  Schmidt  number).  The  flow  term  (K1  Re  Sc  '  ) 
is  dominant  over  all  the  pan  sizes  considered  here.  However,  by  going 
to  very  small  pans  (-  1/2  cm.  rad.  or  less)  the  pure  diffusion  term 
woul  dominate.  The  average  evaporation  rate  per  unit  area  increases 
tremendously  as  one  goes  to  very  small  pans  and  the  average  velocity  of 
the  vapor  leaving  the  liquid  surface  eventually  becomes  large  compared 
to  the  background  convective  air  speed  in  the  room.  Terms  accounting 
for  vapor  density  variations  have  not  been  included  and  do  not  appear 
to  be  important  for  predicting  the  evaporation  rate  under  the  conditions 
considered  here. 

For  moderate  values  of  the  convective  air  speed  the  combined  rate 

3/2 

expression  pr,  :icts  a  nearly  a  dependence  of  the  evaporation  rate 

l 

for  a  particular  liquid  (—  =  0) .  This  is  in  close  agreement  with  ex- 

cl 

perimental  results  obtained  years  ago  by  Ferguson  and  Thomas14  for  water 
evaporating  into  a  laboratory  room  (S./'a  «  0)  . 
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VII.  BINARY  MIXTURES  OF  ORGANIC  LIQUIDS 


The  mass  evaporation  rate  for  a  mixture  of  two  liquids  can  be 
written  as  follows 
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relatively  small  (<  -  .2) .  The  evaporation  rate  is  essentially  a  sum 
of  the  rates  of  the  individual  constituents  of  the  mix  with  provision 
made  for  the  non -ideal  behavior  of  the  mix  through  inclusion  of  the 
activity  coefficients  (y's).  The  x's  are  the  liquid  mole  fractions. 
The  liquid  is  assumed  to  be  thoroughly  mixed. 


2* 


Because  of  the  possible  non-ideal  behavior  of  the  mix  (Raoult's 
law  not  obeyed)  it  is  possible  for  a  mixture  of  two  liquids  to  evaporate 
faster  than  either  of  its  two  constituents.  This  fact  can  be  demonstrated 
using  the  above  equation  and  the  fact  that  the  activity  coefficients 
may  be  greater  than  unity  for  certain  mixtures .  Such  a  non-ideal 
mixture  is  formed  by  cyclohexane  and  ethanol,  and  experimental  results 
demonstrating  that  the  mixture  evaporates  more  rapidly  than  cyclohexane 
(che  more  rapidly  evaporating  component)  are  shown  in  Figure  9. 

A  computer  program  has  been  written  for  the  HP  9810A  computer  to 
calculate  mass  evaporation  rate  as  a  function  of  time  for  binary 
mixtures.  Other  quantities  of  interest  such  as  average  liquid  molecular 
weight,  fc/a,  liquid  mole  fractions,  and  liquid  volume  can  be  printed 
out  as  a  function  of  time  also.  Quantities  assumed  constant  are  the 
liquid  surface  temperature,  viscosity,  and  the  average  diffusion 
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Figure  9.  Mass  Evaporation  Rate  of  the  Binary 
Azeotropic  Mixture,  Cyclohexane  plus 
Ethrmol . 


coefficient.  These  are  not  bad  assumptions  for  the  cases  considered  so 
far.  The  liquid  is  assumed  to  be  well  stirred  which  is  probably  a  good 
assumption  for  thin  liquid  layers.  The  code  predicts  that  thin  layers 
of  mixtures  lose  their  volatiles  faster  than  thicker  layers. 

A  mixture  of  one-third  pentane  and  two-thirds  octane  by  volume  was 
considered  as  a  trial  case  for  the  program.  Twelve  milliliters  of 
liquid  were  placed  in  a  pan  3.7  cm  in  radius  and  0.57  cm  deep.  The 
experimental  results  are  shown  in  Figure  10  together  with  the  theoretical 
curve  for  comparison.  Note  that  the  theory  gives  excellent  agreement 
over  a  range  from  one  to  two  minutes  up  to  ten  minutes  or  perhaps  longer. 
In  the  early  stages  the  experimental  data  is  thought  to  be  higher  be¬ 
cause  the  liquid  surface  temperature  is  slightly  higher  than  the  median 
value  stated.  Also  a  short  time  is  required  to  achieve  pseudo-steady 
state  vapor  concentration  profiles  leading  to  a  higher  initial  rate. 

Binary  mixtures  can  be  used  to  approximate  the  vaporization 
characteristics  of  common  vehicle  fuels.  Figure  10  shows  how  closely 
the  one-third  pentane  two-thirds  octane  mixture  approximates  a  sample 
of  combat  gasoline  vaporizing  from  an  open  pan.  Agreement  is  seen  to 
be  quite  good  from  three  to  ten  minutes  and  probably  longer.  Addition 
of  another  less  volatile  component  would  be  expected  to  extend  the 
agreement  to  very  long  times.  Addition  of  a  more  volatile  component 
would  be  expected  to  improve  agreement  at  short  times . 

VIII.  FLAMMABLE  SURFACES  DEVELOPMENT  FOR  PURE  LIQUIDS, 

MIXTURES  AND  FUELS 

The  time  development  and  extent  of  the  flammable  volume  over  a  fuel 
spill  are  important  input  for  the  military  fuel  fire  problem. 

The  steady  state  flammable  volume  over  a  vaporizing  pure  hydrocarbon 
(hexane)  in  a  circular  pan  is  shown  in  Figure  7  [£/a  u  0.08],  The  experi 
mentally  determined  surface  shown  in  the  figure  represents  the  lower 
flammable  limit  (LFL)  concentration  of  vapor  in  air.  The  upper  flammable 
limit  (UFL)  concentration  is  much  closer  to  the  liquid  surface.  The 
b'-”*4-  '  xtent  of  the  flammable  volume  depends  critically  on  i/a  for 
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Figure  10.  Evaporation  Rates,  Combat  Gasoline  and 
1/3  Pentane -2/3  Octane  Mixture. 


small  l/a.  A  large  box  surrounds  the  pan  and  limits  the  influence  of 
external  air  currents  on  the  system.  The  air  motion  in  the  box  is 
random  with  an  average  air  speed  of  some  2-3  cm/sec  (DISA  and  TSI  low 
velocity  anemometers  employed,  velocity  slightly  lower  than  in  the 
laboratory  room) . 

Slight  changes  in  time,  of  the  velocity  distribution  representing 
the  air  motion  in  the  vicinity  of  the  pan,  imply  that  the  flammable 
envelope  will  fluctuate  with  time.  The  average  flow  direction  and 
magnitude  vary  slightly  with  time  to  produce  small  changes  in  the  extent 
of  the  envelope.  Measurements  show  that  the  hexane  envelope  has  a 
radial  uncertainty  in  the  plane  of  the  pan  of  ±  1  cm  for  vaporization 
from  the  3.7  cm  pan  under  the  flow  conditions  in  the  box.  As  the  flow 
velocity  is  increased  in  one  direction  parallel  to  the  liquid  surface 
the  flammable  envelope  becomes  increasingly  elliptical  in  shape  as  viewed 
from  above.  A  capacitance  spark  probe  is  used  to  locate  the  flammable 
limit  surfaces.  The  entire  apparatus  is  shown  in  Figure  11.  The  spark 
energy  (several  hundred  mj)  is  much  above  the  minimum  ignition  energy 
for  the  particular  electrode  spacing  (~  l  mm)  employed,  according  to  the 
criteria  of  Lewis  and  von  Elbe  as  discussed  in  Mullins  and  Penner15. 

The  experimentally  determined  vertical  height  of  the  LFL  surface 
for  hexane  over  a  3.7  cm  radius  pan  is  approximately  1.0  cm  (for  1/ a 
a  0) .  Three-dimensional  diffusion  theory  (Section  V)  would  predict  a 
much  higher  value,  on  the  order  of  30  cm.,  showing  the  inadequacy  of  a 
pure  diffusion  approach  at  this  relatively  large  pan  size.  However,  for 
smaller  pan  size  pure  diffusion  effects  become  more  important  due  to 
increasing  per  unit  area  vapor  flux  at  the  liquid  surface  (see  combined 
rate  expression  for  small  Re).  For  the  3.7  cm  and  larger  pans  an 
approximate  rule  has  been  developed  for  predicting  the  maximum  height 
of  the  LFL  curface.  For  these  pan  sizes  the  flow  velocity  in  the  box 
is  dominant  in  determining  the  vertical  height  of  the  LFL  surface.  It 
has  been  found  that  the  vertical  height  predicted  from  the  2D  boundary 
layer  flow  theory  at  a  distance  'a*  (one  pan  radius)  from  the  leading 
edge  agrees  approximately  with  the  measured  heights  at  the  center  of 
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the  3.7  cm  and  larger  radii  pans.  For  example,  for  the  3.7  cm  pan  the 
computed  height  (using  the  LFL  concentration)  is  about  1.1  cm  which 
agrees  with  the  measurement  (assuming  a  3  cm/sec  flow) .  The  2D  theory 
predicts  an  a  '  power  dependence  for  the  vertical  height  which  is 
approximately  what  is  observed  as  one  goes  to  larger  pan  sizes.  This 
same  calculation  appears  to  work  for  the  UFL  surface  also,  but  experi¬ 
mental  data  is  not  complete. 

The  2D  theoretical  results  are  presented  by  Bird,  Stewart,  and 
Lightfoot7.  Computations  of  flammable  heights  based  on  these  results 
depend  on  the  evaluation  of  the  incomplete  gamma  function.  An  interest¬ 
ing  relationship  between  the  incomplete  gamma  function  and  the  confluent 
hypergeometric  function  exists  which  permits  an  accurate  evaluation  of 
the  gamma  function.  The  hypergeoraetric  function  is  tabulated  in  the 
NBS  Handbook  of  Mathematical  Functions.  The  relationship  is  as  follows 


r  (y,  u) 
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)M  <1.  J.  »)  -  H  ul/3 
)  -  )  2  .<>8 


(SI) 


where  r  (~,  u)  is  the  incomplete  r  function  according  to  Bird  and 

J  '  I 

T  (y,  »)  is  the  complete  T  function  (■  2.68).  M  (1,  y,  u)  is  the 
confluent  hypergeometric  function.  The  incomplete  r  function  is  re¬ 
lated  to  the  mole  fraction  vapor  concentration  by 
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where  Xp  is  the  vapor  mole  fraction  at  n,  and  Xj^  is  the  mole  fraction 
at  the  liquid  surface.  Also 
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and 


•y  *  is  the  height  above  the  liquid  surface  at  distance  'x'  from  the 
leading  edge  (»x*  taken  to  be  equal  to  'a'  in  this  case),  u  is  -  3 

CD 

cm/sec,  the  average  air  speed. 

Note  the  great  horizontal  spread  of  the  experimentally  determined 
LFL  surface  (Figure  7)  extending  some  7  cm  beyond  the  pan  rim  (3.7  cm 
radius  pan,  ft/a  *  0.08,  for  hexane).  No  mathematical  theory  is  pre¬ 
sented  to  describe  this  phenomenon  but  it  is  felt  to  be  due  to  the  high 
molecular  weight  of  the  hexane  vapor  which  causes  it  to  fall  close  to 
the  base  plate  after  it  clears  the  pan  rim.  Only  liquids  of  relatively 
high  evaporation  rate  (and  proper  flammable  limit)  exhibit  this  large 
horizontal  spread. 

Thus  far  discussion  has  been  limited  to  the  steady  state  surface. 
Time  development  of  the  LFL  surface  of  hexane  is  shown  in  Figure  12. 

The  time  development  study  was  confined  to  a  study  of  the  horizontal 
development  of  the  surface  (in  the  plane  of  the  pan)  as  this  was  easiest 
to  measure  (ft/a  *  0.08  in  Figure  12).  Note  that  the  flammable  region 
develops  rapidly  at  first  and  then  more  slowly  with  steady  state  reached 
in  some  five  to  ten  minutes  at  -  7  cm.  Two  types  of  initial  conditions 
were  employed.  In  one  the  liquid  was  poured  into  the  pan  (S  seconds 
required  to  pour)  and  time  zero  was  taken  at  the  end  of  the  pour.  In 
the  other  the  liquid  was  poured  into  the  container  then  covered  quickly 
with  a  glass  plate.  The  glass  plate  was  then  removed  to  begin  the  test 
(time  zero).  Both  methods  gave  closely  similar  results. 

Mixtures  and  fuels  showed  interesting  flammable  surface  development. 
Results  (again  in  the  plane  of  the  pan)  are  shown  for  a  1/3  pentane  - 
2/3  octane  mixture  (by  volume)  and  for  a  sample  of  combat  gasoline  in 
Figure  13.  The  results  obtained  for  the  (again  ft/a  «  0.08)  binary 
mixture  and  for  the  gasoline  are  indistinguishable  over  a  period  of 
10  minutes.  Note  the  development,  an  initial  rapid  increase  occurs 
followed  by  a  more  gentle  rise  with  a  maximum  reached  in  some  two  or 
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Figure  12.  Time  Development  of  the  Flammable 
Surface  of  Hexane  (In  the  Plane  of 
the  Liquid  Surface). 
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Figure  13.  Flammable  Surface  Development  for  a 
Typical  Sample  of  Combat  Gasoline 
and  for  a  Mixture  of  1/3  Pentane* 
2/3  Octane. 
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three  minutes  (at  8  or  9  cm.  for  the  LFL  surface) .  Then  the  surface 
starts  to  shrink  as  the  loss  of  volatiles  begins  to  dominate,  h  thin 
layer  of  gasoline  allowed  to  sit  in  an  open  space  for  some  fifteen 
minutes  is  a  good  bit  safer  than  a  freshly  poured  sample. 


IX.  EVAPORATION  RATES  FOR  PURE  LIQUIOS  WITH  UNKNOWN 
SURFACE  TEMPERATURE  (THE  HEAT  TRANSFER  PROBLEM) 

The  theoretical  evaporation  rate  expression  developed  in  Section.  VI 
for  pure  liquids  depended  on  a  knowledge  of  the  surface  temperature  of 
the  evaporating  liquid.  If  the  liquid  surface  temperature  is  unknown 
the  problem  becomes  much  more  complicated  since  the  heat  transfer 
characteristics  of  the  system  must  also  be  evaluated. 


Figure  U  depicts  the  energy  balance  in  pseudo -steady  state  vapor¬ 
isation.  In  the  near  steady  state  the  heat  liberated  by  the  liquid 
surface  (Q^  due  to  evaporation  is  balanced  by  the  heat  supplied  by 
conduction  and  convection  through  the  base,  liquid. and  air  f(-IL)(T  -  T^)) , 
^  is  obtained  by  multiplying  (the  mass  evaporation  flux,  sec  Section 
VI)  by  the  heat  of  vaporisation  L.  is  the  heat  transfer  coefficient 
of  the  air,  liquid  and  base.  In  the  near  steady  state  the  liquid  surface 
will  be  at  some  temperature  (Ts)  less  than  the  ambient  temperature  (T^), 

The  evaporative  beat  flux  curves  far' wo  liquids  are  shown  In 


Figure  14.  The  liquid  with  the  lower  evaporative. heat  flux  4.  (less 

than  4)  has  a  surface  temperature  (f;  )  closer  to  msbient  other 
■l  2 

factors  being  equal,  the  implication  is  that  either  the  evapora¬ 
tion  rate  or  the  heat  of  vaporisation  or  both  are  less  for  liquid  ,W«. 
Simple  qualitative  experiments  easily  demonstrate  these  effects. 


The  heal  transfer  coefficient  (K^)  can  be  altered  by  chafing  the 

base  material  or  viscosity  of  the  liquid  (by  adding  a  thickener!  or 

both.  An  aluminum  or  lead  base  plate  gives  a  high  surface  temperature, 

(T  for  liquid  13*  a  glass  base  gives  a  lower  temperature,  and  a  wood 
*1 


or  paper  base  gives  a  Still  lower  temperature 


for  liquid  1) . 
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figure  14.  Effect  of  Sub-trate  on  Vapori lotion. 


The  temperature  effects  are  most  pronounced  for  liquids  with  high  heat 

flux  rates  (like  0  in  the  figure) .  Increasing  the  viscosity  of  a 
1 

liquid  by  adding  a  thickener  decreases  the  heat  transferred  through  it. 

An  evaporating  liquid  is  in  constant  motion16  and  transfers  a  significant 
amount  of  heat  by  convection.  A  thickener  was  added  to  cyclohexane  to 
restrict  the  motion.  The  result  as  expected  was  a  reduction  in  the 
surface  temperature  with  a  consequent  reduction  in  evaporation  rate 
(on  aluminum  base  plate) . 

X.  EVAPORATION  MEASUREMENTS  WITH  THE  TGA 

The  thermo gravimetric  analyzer  (Dupont  Model  9S0)  is  a  very  useful 
instrument  for  making  rapid  determinations  of  relative  evaporation  rates. 
Fairly  reliable  measurements  can  be  made  in  a  period  of  a  few  minutes 
whereas  other  methods  require  considerably  longer.  The  TGA  test  tube 
section  is  shown  in  Figure  15.  The  test  liquid  whose  evaporation  rate 
is  to  be  measured  is  placed  in  the  boat.  Dry  N2  gas  is  passed  over  the 
liquid  parallel  to  the  surface  (a  glass  ball  flow  meter  is  used  to 
monitor  the  N2  from  a  gas  bottle)  .  The  boat  is  partially  entered  into 
the  furnace,  the  degree  of  entry  is  adjusted  until  the  liquid  surface 
temperature  (as  measured  by  the  thermocouple  which  is  placed  close  to 
the  liquid  surface)  is  steady  at  20°C.  The  temperature  is  monitored  at 
intervals  throughout  the  test.  The  volatility  and  heat  of  vaporization 
of  the  liquids  determine  how  far  the  boat  must  be  pushed  into  the 
furnace.  The  evaporation  rate  is  measured  after  steady  state  tempera¬ 
ture  is  reached  by  taking  the  slope  of  the  mass  loss  time  curve  as 
recorded  :n  the  TGA  chart  recorder. 

The  evaporation  rates  measured  with  the  TGA  correlate  approximately 

with  the  product  of  the  diffusion  coefficient  (D) ,  molecular  weight  (My) , 

and  vapor  pressure  (PpQ)  (Figure  16) .  The  TGA  evaporation  rate  is 

equally  divided  between  pure  diffusion  and  flow  controlled  evaporation. 

2/3 

For  flow  controlled  evaporation  the  correlation  would  involve  D  (in 
place  of  D) ,  such  a  correlation  is  also  good  in  this  case.  Uncertainties 
in  the  TGA  data  are  due  to  uncertainties  in  surface  temperature,  in 


APPROXIMATE  SCALE 

Figure  15.  Thermogravimetric  Analyzer. 
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Figure  16.  Evaporation  Rate  Correlation  Determined  from 
DuPont  950  Thermogrovimetric  Analyzer. 


diffusion  coefficient,  and  in  pan  orientation.  The  change  in  depth  of 
the  liquid  below  the  boat  rim  during  the  course  of  a  run  is  not  suffici¬ 
ent  to  cause  more  than  a  ±  10%  variation  in  the  evaporation  rate.  Evapo 
ration  rate  measurements  were  made  with  the  boat  between  two-thirds  and 
one -third  full. 


XI .  CONCLUSIONS 

Even  very  slow  air  movement  (on  the  order  of  a  few  cm/sec.),  such 
as  might  be  encountered  in  a  laboratory  room,  significantly  affects  the 
vaporization  rate  and  flammable  envelope  of  a  pool  of  liquid  hydrocarbon 
The  effects  of  such  small  air  currents  on  evaporation  from  pools  have 
not  been  generally  recognized  much  less  quantified  up  to  this  time. 

The  mass  evaporation  rate  expression  developed  in  Section  VI  has  been 
shown  to  effectively  model  evaporation  into  slow  moving  air  where  the 
velocity  distribution  is  reproducible  over  the  period  of  the  evaporation 
rate  measurements.  The  expression  successfully  predicts  rates  from 
decane  in  the  smallest  pan  (2  cm.  diameter)  to  acetone  in  the  largest 
pan  (17  cm.  diameter),  a  factor  of  one  thousand  in  evaporation  rate. 

The  evaporation  rate  of  six  other  liquids  intermediate  between  decane 
and  acetone  are  also  correctly  predicted.  The  expression  used  here  has 
general  applicability  for  many  types  of  evaporation  problems. 

Vaporization  rates  of  fuels  are  quite  time  dependent  since  fuels 
are  composed  of  many  different  hydrocarbons.  The  vaporization  charac¬ 
teristics  of  common  vehicle  fuels  have  been  approximated  here  at  BRL 
by  simple  mixtures  of  pure  hydrocarbons.  Even  a  binary  mixture  can  do 
an  effective  job  over  a  moderate  length  of  time  as  has  been  demonstrated 
in  the  laboratory.  A  particular  sample  of  regular  combat  gasoline  has 
been  effectively  modeled  over  a  period  of  at  least  ten  minutes  by  a 
mixture  of  1/3  pentane  -  2/3  octane  (by  volume) .  The  great  advantage 
of  using  a  binary  mixture  is  that  accurate  theoretical  predictions  of 
its  vaporization  rate  (and  hence  that  of  the  fuel)  can  be  made  using  the 
computer  program  described  in  Section  VII.  This  program  could  be  ex¬ 
tended  to  treat  n-component  mixtures  and  hence  do  a  more  effective  job 
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at  approximating  a  fuel. 

The  flammable  region  associated  with  a  spill  of  liquid  fuel  or  any 
combustible  liquid  above  its  flash  point  is  not  limited  to  the  physical 
size  of  the  spill  but  extends  appreciably  beyond  and  above  the  spill 
where  the  vapor-air  mixture  is  flammable.  Air  motion  is  critical  in 
determining  the  extent  of  the  envelope.  Detailed  studies  of  the  time 
development  of  the  flammable  envelope  surrounding  pools  of  various  com¬ 
bustible  liquids  have  been  conducted  here  at  BRL.  For  a  pure  liquid, 
such  as  hexane  exposed  to  slow  moving  air  (see  Figure  1),  the  flammable 
surface  develops  rapidly  at  first  and  then  more  slowly  with  steady 
state  reached  in  approximately  five  minutes  after  the  pool  is  initially 
exposed  to  the  atmosphere.  The  vertical  height  of  the  steady  state 
surface  is  proportional  to  the  square  root  of  the  pan  radius  and  is 
relatively  small  compared  to  the  horizontal  extent  of  the  surface.  For 
a  mixture  of  hydrocarbons  such  as  gasoline  the  flammable  envelope 
expands  at  first, as  with  hexane, but  then  dramatically  collapses  as  the 
more  volatile  components  are  lost.  Thus  a  pool  of  gasoline  exposed  to 
the  open  air  is  relatively  safer  after  a  prolonged  length  of  time  than 
it  is  in  the  first  few  minutes.  The  binary  mixture  1/3  pentane  -  2/3 
octane  shows  a  flammable  envelope  which  develops  in  nearly  identical 
fashion  to  the  sample  of  combat  gasoline.  Knowing  the  appropriate 
binary  mixture  enables  one  to  calculate  the  height  of  the  flammable 
envelope  as  a  function  of  time.  Results  such  as  these  show  how  very 
useful  mixtures  can  be  in  approximating  fuels. 

In  conclusion  then,  advances  have  been  made  in  the  theoretical 
description  of  liquid  pool  evaporation  and  in  the  theoretical  and 
experimental  analysis  of  flammable  envelope  development.  Also  consider 
able  success  has  been  achieved  in  the  use  of  a  simple  mixture  to  model 
a  real  fuel. 
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